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The Lieb-Robinson bound shows the existence of a maximum speed of signal propagation in dis- 
crete quantum mechanical systems with local interactions. This generalizes the concept of relativist ic 
causality beyond field theory, and provides a powerful tool in theoretical condensed matter physics 
and quantum information science. Here, we extend the scope of this seminal result by considering 
general Markovian quantum evolution, where we prove that an equivalent bound holds. In addition, 
we use the generalized bound to demonstrate that correlations in the stationary state of a Markov 
process decay on a length-scale set by the Lieb-Robinson velocity and the system's relaxation time. 
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In relativistic quantum field theory, the actions of an 
observer can only influence his future light-cone. This 
can be seen as a consequence of the fact that interac- 
tions are covariant and local, i.e., they couple the field 
at a given point only to the field at points located in- 
finitesimally close to it. The situation is, at first glance, 
quite different for discrete quantum mechanical systems 
with local interactions, such as spin lattices with nearest 
neighbor couplings. There, it is in principle possible to 
send information between any two connected regions in 
an arbitrarily short time, despite the fact that interac- 
tions are local. However, a result first derived by Lieb 
and Robinson [TO] , and improved in Refs. [4, 5j 8] [13], 
demonstrates the existence of an effective light-cone such 
that the amount of information signaled beyond it decays 
exponentially. 

This result, known as the Lieb-Robinson bound, is de- 
rived under the assumption of unitary evolution, i.e., 
when the dynamics is governed by Schrodinger's equa- 
tion. A more general form of quantum evolution is given 
by Markovian dynamical semigroup equations, that are 
the natural generalizations of stochastic processes to the 
quantum setting. These are needed, for instance, to de- 
scribe systems with dissipation or decoherence, and in- 
clude unitary evolution and classical stochastic evolution 
as special cases. 

In this Letter, we extend the Lieb-Robinson bound to 
general local Markovian dynamics. Moreover, we demon- 
strate that the correlations displayed in the stationary 
state of a Markov process decay exponentially beyond 
a length-scale set by the Lieb-Robinson velocity and the 
system's relaxation time, which can be related to the gap 
of the semi-group generator. This also generalizes the re- 
sults of Hastings [5] established in the setting of classical 
Markovian dynamics. 

There are several motivations to study the existence of 
an effective light-cone under general quantum dynamics. 
First, experimental systems are always subject to some 
amount of dissipation and decoherence, so understanding 
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the origin of causality under these conditions is impor- 
tant. Second, the Lieb-Robinson bound has proven to be 
a powerful tool to characterize the structure of ground 
states of gapped Hamiltonians. For instance, Hastings 
and collaborators have used this bound to rigorously 
prove the stability of topological order [2 , the existence 
of PEPS representation of ground states [6], the exponen- 
tial decay of correlations [5], and to generalize the Lieb- 
Schultz-Mattis theorem to higher dimensions [4 . One 
might naturally expect similar characterizations of ther- 
mal states to emerge from the current work; in fact, we 
provide a first step by proving clustering of correlations 
for the fixed points of gapped Markov processes. 

Thirdly, the Lieb-Robinson bound for unitary pro- 
cesses is an important tool in quantum complexity theory, 
see e.g. [7]. It was recently demonstrated that dissipation 
is a universal resource for quantum computation [16]; our 
results complement this finding in a natural way. One 
major open question in this field, related to the quantum 
PCP conjecture, is to identify the complexity of finding 
the ground state energy density of a local Hamiltonian 
within constant accuracy. A problem that is at least as 
hard can be formulated in terms of thermal states, so our 
result could shed new light on this open question. 

Lastly, the existence of a fundamental minimal length 
scale, the Planck length, suggests that physics might 
be fundamentally discrete; many approaches to quan- 
tum gravity have this discreteness built in (e.g. [12] and 
references therein). The black-hole evaporation problem 
also suggests that quantum mechanics could be funda- 
mentally non-unitary, with unitary dynamics emerging 
as a low energy approximation [9] . Our result provides a 
mechanism for emergent causality in such fundamentally 
discrete and non-unitary theories. 

Lieb-Robinson bound — We consider the setting where 
particles are located over a set of vertices A. The parti- 
cle at location x G A has Hilbert space % X1 so the entire 
Hilbert space is ® xe \7~(>x- For any subset of vertices 
X C A, we write Hx — ® x ex^x- We make no dis- 
tinction between operators Ox £ B(Hx) on X and their 
natural embedding in B(7~L). A metric d(x,y) is defined 
between particles locations. A good example to keep 
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in mind are spins located at the vertices of a regular 
D-dimensional lattice, where d(x,y) is the usual graph 
distance. The Hamiltonian is given by a sum of terms 
^ZxcA Hx where Hx = for all X of diameter greater 
than some constant d*. Recall that the diameter of a sub- 
set of vertices is given by the largest distance between any 
pair of vertices inside it. Thus, H is "local" or "short- 
ranged" in the usual sense. 

To motivate the statement of the Lieb-Robinson 
bound, imagine that one observer, Alice, has access to 
some particles A C A and wants to signal to a second ob- 
server, Bob, who has access to B C A. The system is ini- 
tially in the state p G B(H). To send the signal "0", Alice 
does nothing, while if she wants so signal "1" , she applies 
a transformation to the particles in her possession, map- 
ping the state p to p r = p + ie[OA, p] where Oa G B(Ha)- 
To read the signal after some time t, Bob must perform 
a measurement on region B that discriminates between 
the state p(t) and p'(t), where time evolution is gov- 
erned by Schrodinger's equation p(t) = —i[H,p(t)\. If he 
makes a measurement described by the operator G 
B(Hb), the probability that he distinguishes the two sig- 
nals is \Ti{0 B [p(t) - p'(t)]}\ = e\Tr{p[0 B (t), O a ]}\ < 
e l|[Os(t), Oa]\\, where time evolution in the Heisenberg 
picture is governed by 0{t) = i[H, 0(t)]. 

The Lieb-Robinson bound shows that 

\\[0 B (t),0 A ]\\<cV\\0 A \\\\0 B \\ew{- dAB ~ Vt } (1) 

where dAB is the distance between the regions A and £>, 
V = min{|^4|, \B\} is the volume of the smallest of the 
two regions, and c, v, and £ > are constants that de- 
pend only on the microscopic details of the model: the 
interaction strength max^cA ll^x|| 5 the radius of inter- 
actions d*, and the maximal degree of the vertices. Thus, 
signals can only propagate at a finite velocity defining 
an effective light-cone. Outside this cone, the probability 
of detecting a signal falls off exponentially. 

Markov dynamical semigroup equations — We now gener- 
alize the setting by considering a broader class of evolu- 
tion equations. Lindblad has shown [11] that the most 
general differential equation for p that 1) is linear, 2) is 
local in time (Markovian), 3) preserves positivity, and 4) 
preserves the trace must have the form 

p = -i[H, p]+J2 L *P L l ~ \ (LlL a p + P L\L a ) (2) 

a 

where H is a Hamiltonian and L a are any operators. In 
general, the Hamiltonian and the operators L a can be 
time-dependent; our result holds in that case as well but 
we consider time-independent generators for simplicity. 
In the Heisenberg picture, this equation gives 

6 = i[H, Oj+J2 L[OL a - \ {L\L a O + OL\L a ) (3) 

a 

=:C[0). (4) 



It is convenient to adopt a super-operator notation, 
viewing B(H) as a vector space. For O G B(H), we 
use the notation \0)) and denote the Hilbert-Schmidt 
inner product {{0\0')) = TrjO+O'}. Then, we can ex- 
press Eq. p) as \0)) = £|0)), and the formal solution is 
\0(t))) = e*|O(0))). We will often switch between the 
two notations. 

Like in the original setting, we are interested in the case 
where the time-evolution generator is given by the sum 
of local pieces, C = ^2 X ca ^-x w ^h ^x = for all X of 
diameter greater than some constant d* . Each term Cx G 
B(B(Hx)) m that sum has the Lindblad form Eq. ^ 
with Hx,Lx, a £ B{T~Lx)- We make the assumption flTj 
throughout that ||£x|| < 1 5 which is equivalent to fixing 
the time units. 

Lieb-Robinson bound for Markov processes — We now 
come to the main result, which is a bound on [0#(£), Oa) 
where the dynamics of Ob is governed by a local Markov 
process as described above. Our proof is inspired by 
that of [8j [13] . The main complication comes from the 
fact that quantum dynamical semigroup equations do 
not obey Leibniz rule §- t (0 A B ) + ^fO B + A ^ff, 
and that backward time evolution is norm-increasing. 
Here, we will only present parts of the proof that are 
distinct from the unitary case, other details can be found 
in[8l[T3]. 

We are interested in the quantity f(t) = [Ob^^Oa] 
that we can expressed as f(t) = Ca^IOb)), where Ca 
is the super-operator defined by the action Ca\Q)) '•= 
|[Q, Oa]))- We can write a differential equation for f(t) 

f(t)=C A Ce Ct \0 B )) (5) 
= C A C A e ct \0 B )) +C A C nA e ct \0 B )) (6) 
= C A f(t)+C A C nA e ct \0 B )) (7) 

where we have broken the Lindblad super-operator in two 
parts, C nA = ^2x-.xnA^o and C A = C-C nA , and we 
used the fact that [Ca-, £>a\ — 0- 

It can easily be verified by differentiating that the so- 
lution to this differential equation is 

fit) = e^V(O) + f e c ^-^C A CnAe Cs \0 B ))ds. (8) 

Because e CAt is norm-contracting for t > 0, it follows 
that 

\\f(t)\\ < l|/(0)|| + ||C^|| f \\L nA e Cs \0 B ))\\ds. (9) 
Jo 

We can now recurse. Define the quantity 

„ o( x, t) _ sup E^M (10) 

Teh x || / || 

where hx is the set of super-operators T of the form 
Eq. ([3} with L a and H G B(U X )- With this definition, 
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it follows that for any X H Y = 0, 1) elements of hx and 
Ly commute, and 2) elements of hx annihilate B(Hy)- 
Repeating the steps leading to Eq. [9j we have 

M OB (X,t) <M Os (X,0)+ I M OB (Y,s)ds. 

Y-.Ynx^o ^° 

From here, the arguments of [H [13] can be used to show 
that [Ob(£), Oa] is bounded by Eq. ([I]). Note that, using 
the techniques of [14 , it should be possible to generalize 
this bound to the case that the Lindblad super-operator 
is the sum of local unbounded terms with bounded com- 
mutators. 

Convergence rate — Before we examine the correlations 
generated by quantum Markov processes, a few words 
about their asymptotic properties are in order. The 
Lindblad super-operator can be written in its Jordan 
normal form C = SJS~ X where J = J>O Jd 3 X^j) 
is the Jordan matrix and we choose R(Ao) > R(Ai) > 
. . . > M(Afc). Trace preservation implies that R(Ao) = 0. 
If the Markov process has a unique stationary state tt 
such that (^ix\e Ct = ((7r|, then the first Jordan block is 
one-dimensional, do = 1, and the gap A = — R(Ai) 
is strictly positive [TJ [15]. Asymptotically, the system 
converges to this unique stationary state limt^oo e Ct = 
Sdiag(l, 0, . . . , O)^- 1 = \I))((tt\. The gap A governs the 
rate of convergence to equilibrium: 

e? t -\I))((*\=S©e>' t M dj S- 1 (11) 

3>0 

where Md is the d x d matrix with l's on the diagonal 
and on its fcth upper diagonal. Because || < e, we 
conclude that \\e Ct - \I))((tt\\\ < ||<S|| 2 e- A * +1 (note that 
we can always choose S such that ||5|| = ||<S _1 ||, which 
we assume henceforth). Thus, the inverse gap of C sets 
the relaxation rate, but the pre-factor ||<S|| 2 can scale 
with the system size in the case that S is ill-conditioned. 
Finding conditions that make the conditioning number 
of S constant for local C is an interesting question that 
we leave open. 

Clustering of correlations — We now demonstrate that, 
when the system has a relaxation time r that is indepen- 
dent of the system size, the fixed state tt exhibits cluster- 
ing of correlations in the sense that (OaOb) ~ (Oa)(Ob) 
for operators supported on regions dAB ^> vt apart. This 
will occur for instance when C is gapped and ||<S|| con- 
stant. 

Starting in any initial state p, the system reaches the 
stationary state tt in time t ~ A -1 , so 

Tv{[7r-p(f)]0 A B } < ||0 A ||||0 B ||||5|| 2 e - A *. (12) 

In particular, we can choose p to be a product state, 
i.e one without any correlations at all. To gain some 
intuition, we shift to the Heisenberg picture, where we 
know that both Oa{P) an d Osit) grow in space at a 
speed v. Hence, provided that regions A and B are 




FIG. 1: The region R is the union of two membranes sur- 
rounding regions A and B respectively. If we turn off C on 
R, the regions A and B are dynamically decoupled. 



separated by dAB ^ vA _1 , the operators Oa{^) and 
Os(t) will still be supported on disjoint regions by the 
time the system equilibrates, so Tr{pOA(t)OB(t)} ~ 
Tr{pO^(£)}Tr{pO#(£)} for any product state p. 

The problem with this intuitive argument is that, due 
to the failure of Leibniz' rule, the operator (OaOb)^) — 
solution to the differential equation j^X{t) = C[X] with 
X(0) = OaOb — is not equal to Oa(£)Ob(£). The crucial 
observation however is that Leibniz' rule holds for any op- 
erators Oa{P) an d Os(t) contained on regions separated 
by at least the interaction range d* . In that case, we 
have C[0 A {t)0 B (t)) = C[0 A {t)]0 B {t) + A (t)£[0 B (t)}. 
Our generalized Lieb-Robinson bound shows that O^(t) 
and Os{t) remain inside their respective light-cones, save 
for an exponentially decaying tail, so the approximation 
(OaOb)^) ~ OA(t)OB(i) is valid for short times. 

Rigorously, consider the region R that is the union of 
two membranes of thickness 2<i*, the first surrounding 
region A at a distance from A, and the second 

surrounding region B in a similar manner, see Figure 
[l] We write the Lindblad super-operator as the sum of 
two terms, the part supported on i?, Cr = ^2 X cR^ x > 
and the rest C n ^ = C — Cr. Define C(rf) — C n R + 
tjCr, such that C(l) = C and C(0) is the Lindblad super- 
operator obtained by turning off all terms supported on 
R. It is clear that the evolution generated by £(0) cannot 
correlate regions A and B because any operators on those 
regions remain confined inside the regions enclosed by the 
membrane. Thus, for any initial state p in which the two 
regions enclosed by the membranes are not correlated — 
such as a product state — we have ((p|e^ ^|OAO#)) = 
ip^^OA^W^OB)) for all A e B{U A ), B e 
B(H B ), and all t. 

Using the integral representation 

e ct = e c(0)t + f 1 I* e ^r,)(t-0)c R e c ^d/3d V , (13) 
J o Jo 

we can express the time-t correlation ((p\e Ct \OAOB)) as 
the sum of two terms. The first ((ple^^O^O^)) dis- 
plays no correlations as explained above. The second can 
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be bounded using the generalized Lieb-Robinson bound: 

(14) 



< 



Jo Jo 
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< I cT/||£ fl ||||0 A ||||0 B ||exp 

10 

< C y||£fl||||OA||||OB||-exp 



drjdp (15) 
d A B - Z vfi 



2e 



^d/3 (16) 
} (17) 



where V = min{|A| + |i?|} and the other constants 
are as in Eq. (fil. For a D-dimensional regular lattice, 
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< cd 



D-f 

AB J 



where c is a constant that depends on 



the microscopic details of the model. In general, we will 
find = poly (dAB)- Combining this bound with 

Eq. (l2\ yields the desired result 



(O a Ob)-(Oa)(Ob) 
= \\O a \\\\O b \\o(\\S\\ 
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V\\£ f 
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e 2 ^ 



(18) 
") (19) 

(20) 



with ii vA- 1 + f. 

We note that some fairly loose bounds have been used 
in this derivation and a tighter bound may be achievable. 
In particular, we ignored the fact that some initial states 
p reach equilibrium much more rapidly than others. We 
could optimize this choice to improve the bound Eq. ( 12 ), 



subject to the constraint that the two regions enclosed by 
the membranes be initially uncorrelated. A natural guess 
would be to choose the tensor product of the marginals of 



7r over the three regions delimited by the membrane. This 
choice could perhaps compensate for an ill-conditioned 
<S. We note however that there appears to exist some 
local gapped Markov model with long-range correlations 
[3]. Thus, the dependence of the correlations on ||5|| may 
be unavoidable. 

Conclusion — The principle of relativistic causality is a 
pillar of modern physics. The Lieb-Robinson bound 
shows that the principle extends beyond relativistic 
quantum field theory, to the setting of discrete quan- 
tum systems with a Hamiltonian that is the sum of lo- 
cal pieces. Here, we have generalized the Lieb-Robinson 
bound by considering a broader family of dynamical sys- 
tems, namely local Markovian quantum evolution, that 
include unitary quantum evolution and classical stochas- 
tic evolution as special cases. The proof of the Lieb- 
Robinson bound under these conditions differs from the 
original one due to the breakdown of Leibniz' rule and of 
the group properties of the time-evolution operator. 

We have used our generalized bound to demonstrate 
that the correlations displayed in the fixed point of a 
Markov process decay exponentially on a length-scale set 
by the system's equilibration time and the Lieb-Robinson 
velocity. While the latter depends only on the micro- 
scopic details of the model, the former can in general 
scale with the system size even when the generator of 
the Markov process is gapped. Describing conditions un- 
der which the equilibration time is set by the gap of the 
generator remains an interesting open question. 
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